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INTRODUCTION 

A gray tone picture possesses some ambiguity within the pixels due to the 
possible multivalued levels of brightness. The incertitude in an image may arise 
from grayness ambiguity or spatial (geometrical) ambiguity or both. Gray"®? 8 
ambiguity means "indefiniteness" in deciding a pixel as white or black. Spatial 
ambiguity refers to "indefiniteness" in shape and geometry of a region e.g., where is 

the boundary or edge of a region? or is this contour "sharp”? 

When the regions in a image are ill-defined (fuzzy), it is natural and also 
appropriate to avoid committing ourselves to a specific (hard) decision e.g., 
segmentation/thresholding and skeletonization by allowing the segments or skeletons 
or contours, to be fuzzy subsets of the image. Similarly, for describing and 
interpreting ill-defined structural information in a pattern (when the pattern m- 
determinary is due to inherent vagueness rather than randomness), it is natural to 
define primitives and relation among them using labels of fuzzy set 

The present article provides various uncertainty measures arising from grayness 
ambiguity and spatial ambiguity in an image, and their possible applications as 
image information measures. The first part of the article consists of definitions of an 
image in the light of fuzzy set theory, and of information measures (arising from 
fuzziness) and tools relevant for processing/ analysis e.g., fuzzy geometrical 
properties, correlation, bound functions and entropy measures. The second part 
provides formulation of algorithms along with management of uncertainties 
(ambiguities) for segmentation and object extraction, and edge detection. The output 
obtained here is both fuzzy and nonfuzzy. Ambiguity in evaluation and assessment of 
membership function has also been described here. 


* Dr Pal is on leave from the post of Professor in the Electronics and 

Communication Sciences Unit, Indian Statistical Institute, Calcutta 700035, 

India 
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IMAGE DEFINITION 

An image X of size MxN and L levels can be considered as an array of fuzzy 
singletons, each having a value of membership denoting its degree of brightness 
relative to some brightness level 1,1 = 0, 1, 2, ... L - 1. In the notation of fuzzy 
sets, we may therefore write 

^ = {M-x( x rrm) = Hmn/ X mni m = 1, 2 . , . M; n = 1, 2, . . . Nj ( 1 ) 
where Hx( x mn)or Pmn/Xmn, (0 < £ l) 

denotes the grade of possessing some property (e.g., brightness, edginess, 

smoothness) by the (m,n)th pixel intensity x^. In other words, a fuzzy subset of 

an image X is a mapping p from X into [0, 1], For any point p e X , p(p) is 
called the degree of membership of p in p. 

One may use either global or local information of an image in defining a 
membership function characterizing some property* For example, brightness or 
darkness property can be defined only in terms of gray value of a pixel Xj^ whereas, 
edginess, darkness or textural property need the neighborhood information of a pixel 
to define their membership functions. Similarly, positional or co-ordinate 
information is necessary, in addition to gray level and neighborhood information to 
characterize a dynamic property of an image. 

Again, the aforesaid information can be used in a number of ways (in their 
various functional forms), depending on individuals opinion and/or the problem to 
his hand, to define a requisite membership function for an image property. 

MEASURES OF FUZZINESS AND IMAGE 
INFORMATION 


The definitions of various entropy and other related measures which represent 
grayness ambiguity in an image (based on individual pixel as well as a collection of 
pixels) are listed below. 


Linear Index of Fuzziness 

Tl (X) = (2/MN)III Mfnn -n nm l 

m n ^ 

= (2/MN)I I minfUron, l-H™) 

N 


m = 1, 2, . . . M; n = 1, 2, . 
Quadratic Index of Fuzziness 


Yq(X) = (2/VMN) 


M^mn } 


n0.5 


m n 


( 2 ) 


(3) 
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Entropy 

H(X) = (l/MNln2)IISn(p mn ) 
m n 

with Sn(M-mn) = ' M-mn In Mmn 0 Mmn ) ln(l * Mmn) 

m = 1, 2, . . . M; n = 1, 2, ... N 

a denotes the degree of possessing some property p by the (m, n)th pixel 
r mn 

x mn . M'mn denotes the nearest two tone version of ]i 

r* mn 


rth Order Entropy 

H r (X) = (- l/k)X |ix(s7 ) log{p(s? )}+{l- p(s? )} logjl - p(s?)}} (5) 

i = 1, 2, ... k 


s[ denotes the ith combination (sequence) of r pixels in X. k is the number of such 

sequences. p(s[) denotes the degree to which the combination si, as a whole, 
possesses the property p. 

Hybrid Entropy 

H hy (X) = - P w log E w - Pb log E b (6) 

. . E w = (l/MN) XX Pmn ex P(l — M-mn) 
with mn 

E b =(VMN)XX(l-p mn )exp{p 

mn) 

m n 

m = 1, 2, . . . M; n = 1, 2, . . . N 
p denotes the degree of "whiteness" of (m, n)th pixel. P and P. denote 

inn 

probability of occurrences of white (p^ 1) and black (p^ 0) pixels respectively. 
E w and E b denote the average likeliness (possibility) of interpreting a pixel as white 
and black respectively. 


Correlation 


C(pi, P2) = 1 ~ 4 


X X {Plmn — M2mn ) 
_m n 



with 

and 


C(pi,P2)=l ifXj + X 2 = 0 

Xi= XX{2pinm-l} 2 
m n 

X 2 = XX{2p2mn-l} 2 

m n m = 


1, 2, . . . M; n = 1, 2, ... N 


( 7 ) 
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^2 ) denotes the correlation between two properties jij and (defined 
over the same domain). and ^ enote degree of possessing the 
properties p.j and respectively by the (m, n)th pixel. 

These expressions (equations 2-7) are the versions extended to two dimensional 
image plane from those defined for a fuzzy set For example, index of fuzziness was 
defined by Kaufmann [1], entropy by DeLuca and Termini [2], rth order entropy and 
hybrid entropy by Pal and Pal [3], and correlation by Murthy, Pal and Dutta 
Majumdar [4]. 

Interpretation 

Let us describe the properties of these measures along with their relevance to 
image processing/ analysis problems. Index of fuzziness reflects the ambiguity 
present in an image by measuring the distance between its fuzzy property plane and 
the nearest ordinary plane. The term "entropy", on the other hand, uses Shannon's 
function in the property plane but its meaning is quite different from the one of 

classical entropy because no probabilistic concept is needed to define it rf(X) gives a 
measure of the average amount of difficulty in taking a decision on any subset of size 

r with respect to an image property. If r = 1 , H r (X) reduces to (unnormalized) H(X) 
of equation (4). Equation (5) is formulated based on the logarithmic behavior of gain 
function. Similar expression using exponential gain function [14] can also be 

defined. Hhy (X) represents an amount of difficulty in deciding whether a pixel 

possesses certain properties or not by making a prevision on its probability of 
occurrence. In absence of fuzziness (i.e.,with proper defuzzification), reduces to 

two state classical entropy of Shannon, the states being black and white. Since a 
fuzzy set is a generalized version of an ordinary set, the entropy of a fuzzy set 
deserves to be a generalized version of classical entropy by taking into account not 
only the fuzziness of the set but also the underlying probability structure. In that 

respect, can be regarded as a generalized entropy such that classical entropy 

becomes its special case when fuzziness is properly removed. 

All these terms, which give an idea of 'indefiniteness* or fuzziness of an image 
may be regarded as the measures of average intrinsic information which is received 
when (Hie has to make a decision (as in pattern analysis) in order to classify the 
ensembles of patterns described by a fuzzy set. 

T(X) and H(X) are normalized in the interval [0, 1] such that 


Pt 1* Ymin — Hmin ** 0 M-mn — 0 for all (m,n)(X) (8a) 

Et 2: y max ~ Umax = 1 for Jijnn = ^.5 for all (m,n) (8b) 

Pr3: y(X) £ Y(X*)(or, H(X)>H(X*)) (8c ) 

arrf Pr4: y(X) = y(x)(or, H(X)>H(x)) (&J) 


where X* is the 'sharpened' or 'intensified' version of X such that 
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Mx *( x mn) — M-x( x mn) ^ M’x( x mn) — 0-5 

and Hx*( x mn )* M'X ( x mn ) if M x nm)-0-5 ( 9 ) 

In other words, Y(x) or H(X) increases monotonically with p, reaches a 
maximum at p = 0.5 and then decreases monotonically. This is explained in Fig. 1. 
It is to be mentioned here that the definition of fuzzy entropy [14,16] based on 
exponential gain function also satisfies the aforesaid properties. 



Figure 1 Variation of Fuzziness with p. 

According to property 8(c), these parameters decrease with contrast enhancement 
of an image. Now through processing, if we can partially remove the uncertainty on 
the grey levels of X, we say that we have obtained an average amount of information 

given by 5y = y(X)- y(X *) or 6H = H(X)-H(X*) by taking a decision bright or 
dark on the pixels of X. The criteria y(X*)< y(X) and H(X*)< H(X) , in order to 
have positive 8y and 5H -values, follow from Eq. (8c). If the uncertainty is 

completely removed, then y(X *) = H(X *) = 0- In other words, Y(X) and H(X) can 
be regarded as measures of the average amount of information (about the grey levels 
of pixels) which has been lost by transforming the classical pattern (two-tone) into a 
fuzzy pattern X. 

It is to be noted that Y(x) or H(X) reduces to zero as long as p^ is made 0 or 1 

for all (m, n), no matter whether the resulting defuzzification (or transforming 
process) is correct or not. In the following discussion it will be clear how H hy , 

takes care of this situation. 

H r (X) has the following properties: 

Pr 1: H r attains a maximum if p^ = 0.5 for all i. 

Pr 2: H r attains a minimum if Pj = 0 or 1 for all i. 

Pr 3: H r > H* r , where H* r is the rth order entropy of a sharpened version 
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of the fuzzy set. 

Pr 4: H r is, in general, not equal to H r , where H r is the rth order entropy 
of the complement set. 

Pj. 5 . H r <, H r+1 when all pi e [0.5, 1] 

H r > H r+1 when all pj e [0, 0.5]. 

Note that the property P4 of equation 8(d) is not, in general, valid here. The 
additional property Pr 5 implies that H r is a monotonically nonincreasing function 
of r for p j e [0, 0. 5] and a monotonically nondecreasing function of r for 

Pi € [0.5, 1] (when 'min' operator has been used to get the group membership 
value). 

When all the Pj values are same, H^X) = H 2 (X) = . . . = H r (X). This is 

because of the fact that the difficulty in taking a decision regarding possession of a 
property on an individual is same as that of a group selected therefrom. The value of 

H r would, of course, be dependent on the p. values. 

Again, the higher the similarity among singletons the quicker is the convergence 
to the limiting value of H r . Based on this observation, let us define an index of 
similarity of supports of a fuzzy set as S = H^/H^ (when = 0, is also zero and 
$ is taken as 1). Obviously, when p, e [0.5, 1] and the min operator is used to 
assign the degree of possession of the property by a collection of supports, S will lie 
in [0, 1] as H r < H r+1 . Similarly, when pj e [0, 0.5] S may be defined as H 2 /H 1 
so that S lies in [0, 1], Higher the value of S the more alike (similar) are the 
supports of the fuzzy set with respect to the property P. This index of similarity can 
therefore be regarded as a measure of the degree to which the members of a fuzzy set 
are alike. 

Therefore, the value of conventional fuzzy entropy (H 1 or Eq. 4) can only 
indicate whether the fuzziness in a set is low or high. In addition to this, the value 
of HT also enables one to infer whether the fuzzy set contains similar supports (or 
elements) or not. The similarity index thus defined can be successfully used for 
measuring interclass and intraclass ambiguity (i.e., class homogeneity and contrast) 
in pattern recognition and image processing problems. 

The aforesaid features are explained in Table 1 when pj e [0.5, 1], min operator 
is used to compute group membership and k in Eq. 5 is considered to be 
10 Cr ,r = l, 2 6. 

Hhy (X) has the following properties. In the absence of fuzziness when MNP b 
pixels become completely black (Pmn = 0) and MNP w pixels become completely 
white (Pmn = l) then E w = P , E^= P^ and Hj, y boils down to two state classical 
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entropy 

H c = — Pw log P w — log I^j , (10) 

the states being black and white. Thus, Hhy reduces to H c only when a proper 


Table 1 : Higher Order Entropy 


lease p x 

HI 

H2 

H3 

H4 

H5 

H6 

s 

1 

(l.i. 1,1,1, l.i, 1,1.1} 

0 

0 

0 

0 

0 

0 

1 

2 

{.5,.5,.5,.5,.5,.5,.5,.5,.5,.5,.5} 

1 

1 

1 

1 

1 

1 

1 

3 

{1,1, 1,1,1, .5..5,.5,.5,.5} 

.5 

.777 

.916 

.976 

.996 

1 

.642 

4 

{.5, .5, .5,. 5,. 5, .6,. 6, .6,. 6,. 6} 

.980 

.991 

.996 

.999 

.999 

1 

.989 

5 

{.6,.6,.65..9,.9,.9,.9,.9,.9,.915} 

.538 

.678 

.781 

.855 

.905 

.937 

.793 

6 

{.8, .8, .8, .8, .8, .8, .9, .9, .9,-9} 

.538 

.613 

.641 

.649 

.650 

.650 

.878 

7 

{.5,.5,.5,.5,.5,.5,.9,.9,.9,.9J 

.748 

.916 

.979 

.997 

1 

1 

.816 

8 

f.7..7,.7..7,.7,.8..8..8..8,.8} 

.748 

.802 

.830 

.841 

.845 

,846 

.932 


defuzzification process is applied to detect (restore) the pixels. |Hhy H c | can 
therefore be acted as an objective function for enhancement and noise reduction. The 
lower the difference, the lesser is the fuzziness associated with the individual symbol 
and higher will be the accuracy in classifying them as their original value (white or 

black). (This property was lacking with y(X) and H(X) measures (equations 2-4) 
which always reduce to zero irrespective of the defuzzification process). In other 

words, |H h y-H c | represents an amount of information which was lost by 
transforming a two tone image to a gray tone. 

For a given P w and P fc (P w + Pb = 1, 0 < P w ,Pb S l), of all possible defuzzified 

versions, Hj, y is minimum for the one with properly defuzzified. 

If |imn =0.5 for all (m, n) then E w =E fe 

and Hhy = -log(0.5exp0.5) (11) 

i.e., Hhy takes a constant value and becomes independent of P w and P b . This is 
logical in the sense that the machine is unable to take decision on the pixels since all 
it values are 0.5. 

r mn 

Let us consider an example of a digital image in which, say, 70% pixels look 
white, while the remaining 30% look dark. Thus the probability of a white pixel P w 

is 0.7 and that of a dark pixel P fe is 0.3. Suppose, the whiteness of the pixels is not 

constant, i.e., there is a variation (grayness) and similar is the case with the black 
pixels. 

Let us now consider the effect of improper defuzzification on the pattern shown 
in case 1 of the Table 2. Two types of defuzzifications are considered here. In cases 
2-4 all the symbols with p = 0.5 are transformed to zero when some of them were 
actually generated from symbol T. In cases 5-6 of Table 2 some of the p values 
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greater than 0.5 which were generated from symbol 1 (or belong to the white portion 
of the image) are wrongly defuzzified and brought down towards zero (instead of 1). 

In both situations, it is to be noted that |h - Hh y | does not reduce to zero. The case 
7, on the other hand, has all its elements properly defuzzified. As a result, Ej and E Q 

become 0.3 and 0.7 respectively and |Hj, y - H c | reduces to zero. 


Table 2: Effect of wrong defuzzification(with p Q = 0.3 and p t = 0.7) 


C&5G 

E 1 

*9 

H hy 

| H - H hy | 

1 

{.9,.9,.8,.8,.7,.6,.5,.5,.4,.3} 

.620 

.876 

.235 

.375 

2 

{.999,.999,.9,.8,.7,.7,.3,.3,.2,.l} 

.576 

.776 

.342 

.268 

3 

{l,l,l,.99,.9,.9,.l,. 1,0,0} 

.450 

.648 

.542 

.068 

4 

(1,1.1, 1,1, 1.0.0.0.0} 

.400 

.600 

.632 

.021 

5 

{.99,.99,.l,.l,.9,.8,.7,.2,.l,.l} 

.630 

.634 

.456 

.154 

6 

(1,1,0, 0,1, 1,1, 0,0,0} 

.500 

.500 

.693 

.082 

7 

(1,1,1,1,1,1,1,0,0,01 



.611 

0 


There had been some attempts [2, 15] to combine probabilistic entropy and 
possibilistic entropy, they failed to have the aforesaid property of the effect of wrong 
defuzzification. The details of classical entropy measures (e.g., higher order, 
conditional and positional) of an image are available in [14, 16]. 

C(Pj, p^) of equation (7) has the following properties. 

a) If for higher values of Pj(X), (^(X) takes higher values and the converse is 
also true then C(Pj, p^ must be very high. 

b) If with increase of x, both Pj and p^ increase then C(jij, p^) > 0. 

c) If with increase of x, Pj increases and decreases or vice versa then 

CtPj.Pj) < 0. 

d) C(p 1 ,p 1 )= 1 

e) COlj.Pj^COlj.Pj) 

f) C(pj, 1-Pj) = -1 

g) C(p 1 ,p 2 ) = C(p 2 ,p 1 ) 

h) -ISCOij, 

i) COtj.Mj) = -0(1-^, Pj) 

j) C(p 1 ,p 2 ) = C(l-p 1 , I-P 2 ) 
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IMAGE GEOMETRY 

The various geometrical properties of a fuzzy image subset (characterized by 
p.^(x mn ) or simply by p) as defined by Rosenfeld [5,6] and Pal and Ghosh [7] are 

given below with illustration. These provide measures of ambiguity in geometry 
(spatial domain) of an image. 

Area The area of a fuzzy subset p is defined as [5] 

a (^) = -f M- (12) 

where the integration is taken over a region outside which p=4). For p being 
piecewise constant (in case of digital image) the area is 

a(p) = Xp (13) 

where the summation is over a region outside which p=0. Note from equation (13) 
that area is the weighted sum of the regions on which p has constant value weighted 
by these values. 

Example 1 Let p be of the form 
0.2 0.4 0.3 
0.2 0.7 0.6 
0.6 0.5 0.6 

Area a(p) = (0.2+0.4+0.3+0.2+0.7+0.6+0.6+0.5+0.6) = 4.1 


Perimeter If p is piecewise constant, the perimeter of p is defined as [5] 

p(p)= I |p(i)-p(j)|*|A(i, j, k)| 

i, j* k v • 

This is just the weighted sum of the lengths of the arcs A(i, j, k) along which the 
regions having constant p values p(i) and p(j) meet, weighted by the absolute 
difference of these values. In case of an image if we consider the pixels as the 
piecewise constant regions, and the common arc length for adjacent pixels as unity 
then the perimeter of an image is defined by 


p(p)= I |kt(i) - ki( j)| 

>. j 

where p(i) and p(j) are the membership values of two adjacent pixels. 
For the fuzzy subset p of example 1, perimeter is 


(15) 


p(p) = |0. 2 - 0. 4| + 10. 2 - 0. 2| + 10. 4 - 0. 3| + 10.4 - 0.7| 
+10.3 - 0.61 + 10.2 - 0.61 + 10.2 - 0.71 + 10.7 - 0.61 
+|0. 7 - 0. 5| + 10. 6 - 0. 6| + 10. 6 - 0. 5| + 10. 5 - 0. 6| 
=2.3 


Compactness The compactness of a fuzzy set p having an area of a (p) and a 
perimeter of p(p) is defined as [5] 


comp(p) = 


_a(p) 

(pM ) 2 


(16) 
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Physically, compactness means the fraction of maximum area (that can be encircled 
by the perimeter) actually occupied by the object. In non fuzzy case the value of 
compactness is maximum for a circle and is equal to tc/ 4 In case of fuzzy disc, 
where the membership value is only dependent on its distance from the center, this 
compactness value is > n / 4 [6]. Of all possible fuzzy discs compactness is 
therefore minimum for its crisp version. 

For the fuzzy subset p of example 1, comp(p) = 4.1/(2.3*2.3) = 0.775. 


Height and Width The height of a fuzzy set p is defined as [5] 
h(n)= Jmaxurondn 

[17] 

where the integration is taken over a region outside which j x = 0. 

Similarly, the width of the fuzzy set is defined by 
w(n) = JmaxUnmdm 

n (18) 

with the same condition over integration as above. For digital pictures m and n can 
take only discrete values, and since p = 0 outside the bounded region, the max 
operators are taken over a finite set. In this case the definitions take the form 

h(n) = Imaxn mn 

n m KW) 


and 


w(it) = Imaxn mn 
m n 


( 20 ) 


m = 1, 2, . . . M; n = 1, 2, . . . N 

So physically, in case of a digital picture, height is the sum of the maximum 
membership values of each row. Similarly, by width we mean the sum of the 
maximum membership values of each column. 

For the fuzzy subset p of example 1, height is h(p) = 0.4+0.7+0.6 =1.7 and 
width is w(p) = 0.6+0.7+0.6 = 1.9. 


Length and Breadth The length of a fuzzy set p is defined as [7] 

l(n) = max (Jumndn) (2] 

where the integration is taken over the region outside which = 0. In case of a 

digital picture where m and n can take only discrete values the expression takes the 
form 

l(p) = max 
m 

Physically speaking, the length of an image fuzzy subset gives its longest expansion 
in the column direction. If p is crisp, p = 0 or 1; in this case length is the 

maximum number of pixels in a column. Comparing equation (22) with (19) we 
notice that the length is different from height in the sense, the former takes the 
summation of the entries in a column first and then maximizes over different 
columns whereas, the later maximizes the entries in a column and then sums over 
different columns. 


*mn 


( 22 ) 
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The breadth of a fuzzy set p is defined as 

b(p) = max (Jpmndm) (2 3) 

where the integration is taken over the region outside which p = 0. In case of a 
digital picture the expression takes the form 


b(n) = max Xp mn 
n vm 

Physically speaking, the breadth of an image fuzzy subset gives its longest 
expansion in the row direction. If (i is crisp, |jl = 0 or 1; in this case breadth is 

the maximum number of pixels in a row. The difference between width and breadth 
is same as that between height and length. 

For the fuzzy subset ji in example 1, length is l(|i) = 0.4 + 0.7 + 0.5 = 1.6 and 
breadth is b(jx) = 0.6 + 0.5 + 0.6 = 1.7. 



Index of Area Coverage The index of area coverage of a fuzzy set may be 
defined as [7] 

IOAC(p) = ^4^ (25) 

In nonfuzzy case, the IOAC has value of 1 for a rectangle (placed along the axes of 

measurement). For a circle this value is rcr / (2r * 2r) = nf 4 . Physically by 
IOAC of a fuzzy image we mean the fraction (which may be improper also) of the 
maximum area (that can be covered by the length and breadth of the image) actually 
covered by the image. 

For the fuzzy subset |i of example 1, the maximum area that can be covered by 
its length and breadth is 1.6* 1.7 = 2.72 whereas, the actual area is 4.1, so the IOAC 
= 4.1 / 2.72 = 1.51. 

It is to be noted that 1 (X)/h (X) < 1 (26) 

b(X)/w(X)< 1 (27) 

When equality holds for (26) or (27) the object is either vertically or horizontally 
oriented. 

Degree of Adjacency The degree to which two regions S and T of an image are 
adjacent is defined as 


a(S,T)= I 

p g BP(S) 


1 + L(p)-r(q)[ l + d(p) 


(28) 


Here d(p) is the shortest distance between p and q, q is a border pixel (BP) of T and p 
is a border pixel of S. The other symbols are having their same meaning as in the 
previous discussion. 

The degree of adjacency of two regions is maximum (=1) only when they are 
physically adjacent i.e., d(p)=0 and their membership values are also equal i.e., |x(p) = 
r(q). If two regions are physically adjacent then their degree of adjacency is 
determined only by the difference of their membership values. Similarly, if the 
membership values of two regions are equal their degree of adjacency is determined by 
their physical distance only. 
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In the following sections, we will be explaining how the aforesaid measures can 
be used for image segmentation and edge detection problems. 

SEGMENTATION AND OBJECT EXTRACTION 

The problem of grey level thresholding plays an important role in image 
processing and vision problems. For example, in enhancing contrast in a image we 
need to select proper threshold levels from its histogram so that some suitable non- 
linear transformation can highlight a desirable set of pixel intensities compared to 
others. Similarly, in image segmentation one needs proper histogram thresholding 
whose objective is to establish boundaries in order to partition the image spaces into 
meaningful regions. This Section illustrates an algorithm where these various 
information measures can be used to make this selection task automatic so that an 
optimum threshold (or set of thresholds) may be estimated without the need to refer 
directly to the histogram. 

Algorithm 1 

Given an L level image X of dimension MxN with minimum and maximum 

gray vales ^min and / m ax respectively, 

Step 1: Construct the membership plane using the standard S function (equation 
(29)) as 

M-mn = M-W = a, b, c) 

(called bright image plane if the object regions possess higher gray values) 
a Hmn = = 1 " S(£ a, b, c) 

(called dark image plane if the object regions possess lower gray values) 
with cross-over point b and a bandwidth A b. The S function as given below in 
equation (29) is shown graphically in Fig. 2 for an L-level image. 

fmn = ^x( x mn ) = ^( x m n ; a > h, c) = 0, x^ <, a (29a) 

= 2[(x mn ~ a )/(c*~a)] £b (29b) 

= l-2[(x mn -c)/(c-a)] 2 ,bSx mn 5c (29c) 
= 1, Xmn £ C (29d) 

with b = (a+c)/ 2, b-a = c-b = A b. 

The parameter b is the cross-over point, i.e., S(b; a, b ,c) = 0.5. A b is the 
bandwidth. 

Step 2: Compute y(X) ,H(X), Comp(X) and IOAC(X) 

Step 3: Vary b between ^min and £ m ax and select those b for which I(X) 
(where I(X)) denotes one of the aforesaid measures or a combination of them) has 
local minima. Among the local minima let the global one have a cross over point s. 

The level s, therefore, denotes the cross over point of the fuzzy image plane 
U . which has minimum grayness and/or geometrical ambiguity. The plane 
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Step 3: Vary b between ^min and t mAX and select those b for which I(X) 

(where I(X» denotes one of the aforesaid measures or a combination of them) has 

local minima. Among the local minima let the global one have a cross over point s. 

The level s, therefore, denotes the cross over point of the fuzzy image plane 

u. , which has minimum grayness and/or geometrical ambiguity. The p plane 

mn 11111 



0 i b < U1 

^ mn 

Figure 2 Standard S function for an L-level image. 

then can be viewed as a fuzzy segmented version of the image X. For the purpose of 
nonfuzzy segmentation, we can take s as the threshold or boundary for classifying or 
segmenting image into object and background. (For images having multiple regions, 
one would have a set of such optimum p(X) planes). 

The measure I(X) in Step 3 can represent either grayness ambiguity (i.e., y( x ) 
or H(X)) or geometrical ambiguity (i.e„ comp(X) or IOAC(X) or a(S,T)) or both 
(i.e., product of grayness and geometrical ambiguities). 

Faster Method of Computation 


From the algorithm 1 it appears that one needs to scan an L level image L times 
(corresponding to L cross over points of the membership function) for computing the 
parameters for detecting its threshold. Hie time of computation can be reduced 
significantly by scanning it only once for computing its co-occurrence matrix, row 
histogram and column histogram, and by computing p(l), 1 = 1, 2, ... L every time 
with the membership function of a particular cross over point 

The computations of y(X) (or H(X)), a(X), p(X), 1(X) and b(X) can be made 
faster in the following way. Let h(i), i=l,2..L be the number of occurrences of the 
level i, c[ij], i = 1, 2 . . L, j = 1, 2 . . L the co-occurrence matrix and p(i), i = 1, 2.. 
L the membership vector for a fixed cross over point of an L level image X. 

Determine y( x ) , area and perimeter as 


Y(X) = -^ ST(i)h(i) 
MN i=1 

T(i)= min{p(i),l-p(i)} 


(30a) 

(30b) 
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For calculating length and breadth following steps can be used. Compute the 
row histogram R[m, 1], m = 1, . . . M, l = 1 . . L, where Rfm, 1] represents the 
number of occurrences of the gray level 1 in the mth row of the image. Find the 
column histogram C[n, 1], n = 1. . N, 1 = 1 . . L, where C[n, 1] represents the 
number of occurrences of the gray level 1 in the nth column of the image. Calculate 
length and breadth as 

1(X)= max I C[n, 1]. p(l) 
n ] = i 

b(X) = max 2 R[m, 1]. p(l) 

m i=i 

Some Remarks 

The grayness ambiguity measure e.g., y(X) or H(X) basically sharpens the 
histogram of X using its global information only and it detects a single threshold in 
its valley region. Therefore, if the histogram does not have a valley* the above 
measures will not be able to select a threshold for partitioning the histogram. This 

can readily be seen from Equation (30) which shows that the minima of Y(X) 
measure will only correspond to those regions of gray level which has minimum 
occurrences (i.e., valley region). Comp (X) or IOAC(X), on the other hand, uses 
local information to determine the fuzziness in spatial domain of an image. As a 
result, these are expected to result better segmentation by detecting thresholds even in 
the absence of a valley in the histogram. 

Again, com p(X) measure attempts to make a circular approximation of the 
object region for its extraction, whereas, the IOAC(X) goes by the rectangular 
approximation. Their suitability to an image should therefore be guided by this 
criterion. 

Choice of Membership Function 

In the aforesaid algorithm w - 2Ab is the length of the interval which is shifted 
over the entire dynamic range of gray scale. As w decreases, the M-O^) plane would 

have more intensified contrast around the cross-over point resulting in decrease of 
ambiguity in X. As a result, the possibility of detecting some undesirable thresholds 
(spurious minima) increases because of the smaller value of A b. On the other hand, 
increase of w results in a higher value of fuzziness and thus leads towards the 
possibility of losing some of the weak minima. 

The criteria regarding the selection of membership function and the length of 
window (i.e., w) have been reported recently by Murthy and Pal [10] assuming 
continuous function for both histogram and membership function. For a fuzzy set 

"bright image plane", the membership function p: [0, w] — » [0,1] should be such 
that 

i) jx is continuous, p(0) = 0, |i(w) = 1 

ii) p is montominally non-decreasing, and 

iii) |i(x) = 1- jx(w-x) for all x € [0, w] where w>0 is the length of the window. 
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Furthermore, p should satisfy the bound criteria derived based on the correlation 
measure (equation 7). The main properties on which correlation was formulated are 
p : If for higher values of Pj, p 2 takes higher values and for lower values 

of (J.j, p 2 a^ 50 takes lower values then C(Pj, P 2 ) > 0 

P 2 : If M-jT and p^ then C(p r p 2 ) > 0 

P 3 : Ifpjt and^^ then C(Pj,p 2 ) < 0 

[T denotes increases and ^ denotes decreases]. 

It is to be mentioned that P2 and P3 should not be considered in isolation of Pj 

Had this been the case, one can cite several examples when Pj T and p^ but Cfp^ 

p 2 ) < 0 and pj t and p^ but C(p r p 2 ) > 0. Subsequently, the type of 

membership functions which should not be considered in fuzzy set theory are 
categorized with the help of correlation. Bound functions hj and h 2 are accordingly 

derived [1 1]. They are 

hj(x) =0, 0 < x <e (35) 

= x- 6, E< X < 1 

h 2 (x)=x+€, 0 < x < 1- e (36) 

= 1, 1- €< X < 1 

where e = 0.25. The bounds for membership function p, are such that 
hi(x) < p,(x) < h 2 (x) for x e [0,1] . 

For x belonging to any arbitrary interval, the bound functions will be changed 
proportionately. For hi <|i^h 2 , C(hj,h2)^0 ,C(hj,p)^0 and C(h 2 ,p)^0, 



x 

Figure 3 Bound Functions for |i(x). 
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The function p lying in between hj and hj does not have most of its variation 

concentrated (i) in a very small interval, (ii) towards one of the end points of the 
interval under consideration and (iii) towards both the end points of the interval under 
consideration. 

Figure 3 shows such bound functions. It is to be noted that Zadeh's standard S 
function (equation 29) satisfies these bounds. 

It has been shown [10] that for detecting a minimum in the valley region of a 
histogram, the window length w of the p function should be less that the distance 
between two peaks around that valley region. 

H r as an Objective Criterion 

Let us now explain another way of extracting object by minimizing higher order 
fuzzy entropy (equation 5) of both object and background regions. Before explaining 
the algorithm, let us describe the membership function and its selection procedure. 

Let s be an assumed threshold which partitions the image X into two parts 
namely, object and background. Suppose the gray level ranges [1 - s] and {s + 1 - L] 
denote, respectively, the object and background of the image X. An inverse rr-type 
function as shown by the solid line in the Figure 4 is used here to obtain p mn values 

of X. The inverse rc-type function is seen (from Fig. 4) to be generated by taking 
union of S(x ; (s - (L - s)), s, L) and 1 - S(x; 1, s, (s + s - 1)), where S denotes the 
standard S function defined by Zadeh (equation 29). 

The resulting function as shown by the solid line, makes p lie in [0.5,1]. Since 
the ambiguity (difficulty) in deciding a level as a member of the object or the 
background is maximum for the boundary level S, it has been assigned a membership 
value of 0.5 (i.e., cross-over point). Ambiguity decreases (i.e., degree of 
belongingness to either object or background increases) as the gray value moves away 
from s on either side. The p mn thus obtained denotes the degree of belongingness of 

a pixel x mn to either object or background. 

Since s is not necessarily the mid point of the entire gray scale, the membership 
function (solid line if Fig. 4) may not be a symmmetric one. It is further to be noted 
that one may use any linear or nonlinear equation (instead of Zadeh's standard S 
function) to represent the membership function in Fig. 4. Unlike the Algorithm- 1, 
the membership function does not need any parameter selection to control the output. 

Algorithm 2 

Assume a threshold s, 1 < s < L and execute the following steps. 

Step 1 : Apply an inverse n - type function [Fig. 4] to get the fuzzy p mn plane, 

with P mn e [0.5, 1]. (The membership function is in general asymmetric). 

Step 2: Compute the rth order fuzzy entropy of the object Hq and the background 
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Figure 4 Inverse n function (solid line) for computing object and background 

entropy. 

Hg considering only the spatially adjacent sequences of pixels present within the 
object and background respectively. Use the ’min' operator to get the membership 
value of a sequence of pixels. 

Step 3: Compute the total rth order fuzzy entropy of the partitioned image as 

H s = H o + H{j. 

Step 4: Minimize H* with respect to s to get the threshold for object background 
classification. 

7 

Referring back to the Table 1 , we have seen that H reflects the homogeneity 

among the supports in a set, in a better way than H* does. Higher the value of r f the 
stronger is the validity of this fact. Thus, considering the problem of object- 

background classification, H r seems to be more sensitive (as r increases) to the 
selection of appropriate threshold; i.e., the improper selection of the threshold is 

more strongly reflected by H r than H r ^ For example, the thresholds obtained by 

measure has more validity than those by H* (which only takes into account the 
histogram information). Similar arguments hold good for even higher order (r > 2) 
entropy. 

Example 2 

Figures 5 and 6 show the images of Lincoln and blurred chromosome along with 
the histogram. Table 3 shows the thresholds obtained by comp (X) and IOAC (X) 
measures for various window sizes w when Zadeh's S function is used as membership 
function. Lincoln image is of 64x64 with 32 gray levels whereas, chromosome 
image is of 64x64 with 64 gray levels. 


227 




Figure 6(a) Input. Figure 6(b) Histogram. 
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Figure 7(b) Threshold = 32. Figure 7(c) Threshold = 56. 


Table 3 Various Thresholds (* denotes global minimum) 

W 


12 
16 
20 
24 

Threshold produced by measure (Algorithm 2) is 8 for Lincoln image. Some 
typical nonfuzzy thresholded outputs of these images are shown in Figure 7. 
Recently, transitional correlation and within class correlation have been defined [12] 
based on equation (7) for image segmentation which takes both local and global 
information into account Automatic selection of an appropriate enhancement 
operator is available in [17]. 
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IOAC 
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32*46 
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34 
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Comp 
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10 

10 
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12 

10 
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16 
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11 
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EDGINESS MEASURE 

Let us now describe an edginess measure [18,19] based on ^(Equation 5) which 
denotes an amount of difficulty is deciding whether a pixel can be called an edge or 

not. Let be a 3 x 3 neighborhood of a pixel at (x, y) such that 
x,y 

N x,y = K*. y). (x-i.y). (x+i.y). (x,y-i). (*.y+i). (x-i. y-i), 

(x-l. y+l), (x+l,y-l), (x+l,y+l)} (37) 

The edge-entropy, h e of the pixel (x, y), giving a measure of edginess at (x, 
x,y 

y) may be computed as follows. For every pixel (x, y), compute the average, 

maximum and minimum values of gray levels over __ . Let us denote the average, 

x,y 

maximum and minimum values by Avg, Max, Min respectively. Now define the 


following parameters. 

D = max { Max - Avg, Avg - Min} (38) 

B = Avg (39) 

A = B - D (40) 

C = B + D (41) 


A n-type membership function is then used to compute it for all (x, y) 

xy 

3 

e ^x.y > suc h that p(A) = p(C) = 0.5 and p(B) =1. It is to be noted that p xy ^ 0.5. 
Such a p , therefore, gives the degree to which a gray level is close to the average 

x y 

3 

value computed over N Xt y . i n other words, it represents a fuzzy set "pixel intensity 

close to its average value", averaged over N x y , When all pixel values over 

are either equal or close to each other (i.e., they are within the same region), such a 

transformation will make all p = 1 or close to 1. In other words, if there is no 

xy 

edge, pixel values will be close to each other and the p values will be close to 
one(l); thus resulting in a low value of H*. On the other hand, if there is an edge 

(dissimilarity in gray values over N x ) t then the p values will be more away from 

1 1 3 

unity; thus resulting in a high value of H 1 . Therefore, the entropy H over N x y 

£ 

can be viewed as a measure of edginess ( H x ,y ) at the point (x, y). The higher the 

value of H x y , the stronger is the edge intensity and the easier is its detection. As 
mentioned before, there are several ways in which one can define a n-type function 
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1 


0.5 — 


Figure 8 
(solid line) as shown in Fig. 8. 

The proposed entropic measure is less sensitive to noise because of the use of a 
dynamic membership function based on a local neighborhood. The method is also 
not sensitive to the direction of edges. Other edginess measures are available in 
[13,20]. 

CONCLUSIONS 

Various uncertainty and image information measures, as conveyed by entropy and 
fuzzy geometry, have been explained. The problems of object extraction and edge 
detection have been considered, as an example, to demonstrate their applications. 
Uncertainty in membership evaluation in these problems and its management have 
also been explained through bound functions. The measures comp(X) and IOAC(X) 
can also be used for skeleton extraction and medial axis transformation of a gray tone 
image [21]. 
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